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Abstract 

The charged Higgs boson decays — > W A\ and — > W^hi are studied in the framework 
of the next-to Minimal Supersymmetric Standard Model (NMSSM). It is found that the decay rate 
for H± -> W ± A 1 can exceed the rates for the and tb channels both below and above the 
top-bottom threshold. The dominance of — > W^A\ is most readily achieved when A\ has a 
large doublet component and small mass. We also study the production process pp — > H^A\ at 
the LHC followed by the decay — > W^A\ which leads to the signature W^A\A\. We suggest 
that pp — > H^Ai is a promising discovery channel for a light charged Higgs boson in the NMSSM 
with small or moderate tan/3 and dominant decay mode —> W^Ai. This M^ ± j4ij4i signature 
can also arise from the Higgsstrahlung process pp — > W^h\ followed by the decay h\ — > A\A\. It 
is shown that there exist regions of parameter space where these processes can have comparable 
cross sections and we suggest that their respective signals can be distinguished at the LHC by 
using appropriate reconstruction methods. 

PACS numbers: 12.60.Pr, 14.80. Cp 



1 



I. INTRODUCTION 



An attractive extension of the Minimal Supersymmetric Standard Model (MSSM) is the 
Next-to MSSM (NMSSM) in which an additional singlet neutral complex scalar field S is 
added. The presence of this singlet field provides an elegant solution to the \x problem of the 
MSSM. The \i parameter in the MSSM superpotential, which does not break supersymmetry 
(SUSY) and is present when SUSY is unbroken, is completely unrelated to the electroweak 
or SUSY breaking scales. In some models like Supergravity, /i is naturally expected to be of 
the order Mpi anc k. However, the radiative electroweak symmetry breaking conditions require 
the /i parameter to be of the same order as Mz- Such a conflict is called the /i problem 

The superpotential of the NMSSM contains the term \S H u Hd, and the fi term of the 
MSSM which mixes the two doublet fields H u and Hd is not present explicitly. When the 
singlet field acquires a vacuum expectation value < s > of the order of the SUSY breaking 
scale, an effective /i parameter fi e ff = As of the order of the electroweak scale is then 
dynamically generated. Moreover, it has been shown that with the additional singlet Higgs 
field the MSSM fine-tuning (or "little hierarchy problem") problem can be ameliorated in 
regions of the NMSSM parameter space 0, Hj[. 

A charged Higgs boson (H^) appears in any extension of the Standard Model with two 
hypercharge Y=l doublets. Its phenomenology has been extensively studied in both the 
Two Higgs Doublet Model (2HDM) and MSSM. The phenomenology of if ± in the NMSSM 
is similar in many ways to that in the MSSM since no charged singlet fields have been added. 
The increased parameter content of the NMSSM scalar potential compared to that of the 
MSSM permits large mass splittings among the Higgs spectrum, which allows other decay 
modes of H ± to be important which were substantially suppressed in the context of the 
MSSM. In the MSSM the coupling H ± AW (where A is the CP-odd neutral Higgs boson) 
contains no mixing angle suppression but the relation Ma ~ Mh± ensures that the decay 
— > AW is greatly suppressed in most of the parameter space In the NMSSM 

there are two pseudoscalars A 1 and A 2 which are mixtures of the doublet and singlet fields. 
There exists regions in the theoretical parameter space where A\ is predominantly doublet 
and light, and hence the decay H ± — > AiW is unsuppressed. 

The importance of the decay H ± — > A{W in the NMSSM was emphasized in j6[ where 
it was shown that dominance over — > cs,tu is possible and branching ratios close to 
100% can be attained for intermediate values of tan/3. A LHC simulation was performed 
in [?J and concluded that such a decay offers very good detection prospects for H ± if the 
branching ratios of t — > H^b and H — > A{W are sufficiently large. In this work we perform 
a comprehensive scan of the NMSSM parameter space using the publicly available code 
NMHDECAY jsj in order to identify the regions where H ± — > A{W can be sizeable. 

The strength of the coupling H^AAW can also have an application to theproduction of 



via pp — > H ± A\ which has been studied in the CP conserving MSSM |9[,[10| and CP 
violating MSSM [ll]. If the branching ratio for the decay H ± — > A±W were also sizeable 
such a production mechanism would lead a final state of Wbbbb (for Ma 1 > Znib) [lH which 
has been simulated [l2| in the context of the LHC with promising conclusions. This Wbbbb 
signature can also arise from the process pp — > Wh\ — > WA±Ai which was simulated in 



13( and shown to provide a clear signal at the LHC. We compare the magnitude of both 



mechanisms and discuss how they may be distinguished. 

Our work is organized as follows: in section II we present a short review of the Higgs 
sector of the NMSSM; in section III the limits that lead to a light A x in the NMSSM 
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parameter space are listed; in section IV the phenomenology of H is introduced; section 
V contains our numerical results for the branching ratios of — > AiW,h\W and cross- 
sections pp — > H ± A\ — > Wbbbb{W tttt) and pp — > Wh\ — > Wbbbb{W tttt) . Conclusions 
are given in section VI. 



II. A BRIEF REVIEW ON THE HIGGS SECTOR OF THE NMSSM 

For detailed discussions of the Higgs sector of the NMSSM the reader is referred to 



U, [15|, [16|, [13, [18|. In this section we follow the notation of Ref. [19|. The NMSSM Higgs 
sector differs from that of the MSSM by the addition of an extra complex scalar field, S. The 
Higgs fields of the model then consist of the usual two Higgs doublets H u and H d together 
with this extra Higgs singlet. 

In the NMSSM Lagrangian, the extra singlet field is allowed to couple only to the Higgs 
doublets of the model and consequently the couplings of the new field S to gauge bosons 
and fermions will only be manifest via their mixing with the doublet Higgs fields. The 
superpotential of the NMSSM is given by 

W = H/mssm + AS H u H d + i« S 3 . (1) 

where H^mssm is the usual MSSM superpotential and only terms that depend on the singlet 
field are explicitly written. The soft breaking terms for both the doublet and singlet are 
included in V so &'- 

Koft = m 2 H \H u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + [XA x SH u H d + ^kA k S 3 + h.c] , (2) 

The parameters additional to those of the MSSM are: A, k, A\, A k , ms and the vacuum 
expectation value of the singlet field, s, which will generate the effective fj, term given by 
fi e ff = As. As in the MSSM, ms can be fixed by the minimization condition of the scalar 
potential. 

After electroweak symmetry breaking the Higgs spectrum of the NMSSM consists of 
three neutral scalars (hi,h 2 ,h 3 ), two pseudoscalars (Ai,A 2 ) and a pair of charged Higgs 
bosons H ± . In both the CP-odd and CP-even sector the physical eigenstates are ordered as 
M hl < M h2 < M h:j and M Al ^ M A2 . The mass of E ± at tree- level is given by E3], 0: 



M 2 H± =^L(Ax + ks) + M 2 W - A V (3) 
sin 2p 

where tan/5 = v u /v d and v 2 = v 2 +v d . This differs from the corresponding MSSM expression 
in which M A and M#± are strongly correlated and become roughly equal for M A > 140 GeV. 

The CP-odd mass matrix can be obtained as follows: Firstly, as in MSSM one rotates the 
bare fields (3mH u , SmHd, 3mS) into a basis (A, G, QinS) where G is a massless Goldstone 
boson. Then one eliminates the Goldstone mode and the remaining 2x2 CP-odd mass 
matrix in the basis (A, SmS) is given by: 

As 

M P,U = —Ik o( A X + Ks), 

sm p cos p 



M 2 P22 = (2Xk + —)sm2f3v 2 - 3kA k s, 

(2 



M 2 P12 = Xv(A x -2ks). (4) 
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Here A = cos /3SmH u + sin /5SmHd is the CP-odd MSSM Higgs boson while SmS comes 
from the singlet S field. The pseudoscalars fields are further rotated to the diagonal basis 
(Ai, A2) by an orthogonal 2x2 matrix such that: 

A\ = cos 9 aA + sin^Sn^S) 

A 2 = -sin9 1 A + cos9 A ^rn(S) (5) 



where 



M 2 P12 M 2 -M 2 P11 
cos 9 a = — ' , sin 9 a = — 1 (6) 



^M P ^ 2 + (Ml-M% u y ^JM% 12 + (Ml 



M 



2 \2 
P,ll) 



The Higgs boson-gauge boson couplings originate from the covariant derivative of the 
kinetic energy term. Those relevant for our study are described by the following Lagrangian: 

Cvvh,vhh = gm w g VVhi W + »W-h % - gW+ C 9w+ ^~ hi K + ^-A t ) ? H~ + h.c (7) 

where g VVhi = sin (3S a + cos (3S i2 , gw+H-h, = cos f3S a - sin (3S i2 , P u = cos9 A and P 21 = 
— sin^, S and P are orthogonal matrix which diagonalize respectively the CP-even and 
CP-odd scalar mass matrix. ^From the last term in eq. ([7]) one can see that the vertex 
W ± H Zf Ai is directly proportional to Pn i.e. the doublet component of the mass eigenstate 
Ai. Consequently, if A\ is entirely composed of doublet fields this coupling is maximized 
and if A\ is purely singlet the coupling vanishes. 

As in the MSSM one can easily derive the following sum rules: 



^ 9\VWhi 



i=l 



9wwhi + 9w+H-hi + S i3 — 1 i — 1, 2, 3 (8) 

Here is the singlet component of hi. From the second sum rule it follows that if hi is 
purely doublet (S i3 » 0) then the MSSM sum rule, 9\ywhi 9w+H-hi — 1, is recovered 
where hi is entirely composed of doublet fields. Conversely, if hi is purely singlet (Sf 3 ~ 1) 
then one has glvwh- + 9w+H-h ~ ^ anc ^ both hjVV and hiH + W~ must be suppressed, and 
this will present a real challenge for the detection of Higgs bosons. This sum rule will be 
explored in our numerical analysis. 



III. A LIGHT A x IN THE NMSSM PARAMETER SPACE 

The parameter space of the NMSSM can naturally accommodate a light A\ which is of 
great phenomenological interest. To identify such regions it is instructive to examine the 
vanishing limits of the determinant of the mass matrix of the pseudoscalar, which can be 
expressed as: 

3kAs / \ 

DetMl = 2ks 2 A k + 2sA K A x - 3\A\v 2 sin 2/5 . (9) 

sm 2p \ I 

It is then straightforward to identify four distinct cases where DetMp approaches 0: 
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• Case 1: A x -> and A K — > 21 

• Case 2: k -> @, 

• Case 3: A — > 0, 

• Case 4: s — > 0. 

Moreover, it is evident that combinations of these basic cases can also lead to a light A 



The requirement of perturbativity up to the grand unification scale restricts A < 0.8 [23 
Therefore Case 4 (s — > 0) is ruled out since it would lead to a very small fi e ff which is 
excluded by the mass bound for charginos from direct searches. However, if one gives up 
this perturbative requirement up to grand unification scale and considers A 3> 1, as in the 
so-called ASUSY model 24j (which can be realized in the supersymmetric fat Higgs models 
(25|). then Case 4 might be viable. 

The first two limits are related with the discrete symmetries of Higgs potential: one is 



called the R-axion limit with A\ — > and A K — > [21J; the other is called the PQ-axion 
limit with k — > the superpotential eq.fJTJ) and its associated Lagrangian contains an extra 
global 17(1) symmetry [22|. In both cases, these symmetries are spontaneously broken by 
the Higgs vev leading to Pseudo-Goldstone boson in the spectrum. 

At tree-level, in the R-axion limit 2l|, the mass spectra and mixing of the CP-odd Higgs 
sector can be expressed as: 



m, 



l a 2 
tan 9 a 



3s(—kA k sin 2 9a + 



2 sin 2(3 



XA x cos 2 9a) + 0(k 2 A 2 k ,\ 2 AI), 



2\kv 2 
cos 2 9 a 



v sin 2(3 



sm2(3 + 0(K 2 A 2 K ,X 2 Ai), 
0(kA k ,XA x ). 



(10) 



In the R-axion limit scenario, as can be seen from eq. (110p . a light pseudoscalar is obtained 
for small kA k and \A\ or a combination of small kA k and \A\. 
At tree-level, in the PQ-axion limit 22|, one has: 



m A 2 



tan 9, 



3sk(—A k sin 2 9a 



2\A\v 

sin 29 a 

2s 

v sin 2(3 



+ 0(k 2 ) 
+ 0(k). 
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sin 2(3 



As cos 2 9 A ) +0(k 2 



(11) 



It is interesting to see that in eq. (fTTj) the limit k —> gives mA 1 —> 0. This is actually the 
case where the U(l) PQ symmetry is left unbroken in the superpotential. The spontaneous 
breaking of such PQ symmetry by a Higgs vev leads to a massless Goldstone boson, the 
axion. To obtain a light pseudoscalar Ay one needs to introduce a small k which only 
slightly breaks the PQ symmetry. 
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The third case is also related with a discrete symmetry of two Higgs doublet models. In 
this limit one has: 



2As 

vr? M = -3ksA k + 0(A) , 

tanlSU = * {Ax 3 ~a™ )V +0& 2 ). (12) 

When A — > 0, a large value for s is needed to keep /i e // of the order of the electroweak scale. 
In this case A — > 0, and for fi e ff fixed, A\ is mainly doublet and this is the exact MSSM 
limit. 



IV. H± IN THE NMSSM 

In this section we describe the phenomenology of the H ± in the NMSSM and highlight 
its differences with the phenomenology of in the MSSM. The phenomenology of H ± 
in the NMSSM has many similarities with that of H ± in the MSSM (the latter recently 



reviewed in |26|). This is to be expected since the fermionic couplings are identical in the 
two models. The main differences in their phenomenology originate from the possibility of 
large mass splittings among the Higgs bosons in the NMSSM which permits decay channels 
like H ± — > A\W to proceed on-shell p. In the MSSM such a decay can only be open 
for extreme choices of certain SUSY parameters (e.g. for \i > AMsusy (27|) which induce 
large quantum corrections in the effective scalar potential. Moreover, in the NMSSM a light 
CP-even h\ is also allowed and one can have the opening of the decay H ± — > h{W both 
below and above the top-bottom threshold. This latter channel may change the NMSSM 
phenomenological predictions for the charged Higgs with respect to the MSSM [gj. In the 
MSSM the decay — > h\W is also open but the coupling gw+H-hi ~ cos 2 (/3 — a) is 
strongly suppressed when Mh± 3> m/ ll + mw and thus its branching ratio is very small for 
such Mn±. For Mh± < tu^ + my/ and just above the threshold the branching ratio for this 
channel can reach 10% at most for small values of tan/? jij, 0]. 

The phenomenology of H in the NMSSM has received considerably less attention than 
its neutral Higgs sector. In recent years much effort has been focused on establishing a "no- 
lose theorem" at the LHC in which detection of at least one Higgs boson in the NMSSM 



is guaranteed. However, the potential importance of the decay h\ — > A\A\ [2l|,[28| has 
prevented such a theorem being established 19], 29], 3(|, [3l|], Moreover, it has been 
shown that a large branching ratio for h\ — > A\A\ would weaken the LEP bounds for a SM 
like h x in the NMSSM M. 

For Ma ± < [l^ll^ dominance of hi — > A\Ai has the virtue of allowing hi as light as 
90 — > 100 GeV and can realize the "LEP excess scenario" easily. Such values of M/ ll can be 
accommodated in the NMSSM with little fine-tuning, in contrast to the MSSM case where 
considerable fine-tuning is necessary in order to comply with the LEP limit > 114 GeV 
from the Higgsstrahlung channel. However, a large branching ratio for hi — > A\A\ followed 
by AiAi — > 4r is challenging for detection at the Tevatron (see [13]). For the final states 
V2b2r and V4b at the Tevatron, the observation is difficult due to the limited statistics, 



as shown in 35(. At the LHC, by utilizing the central exclusive production process and 



high-resolution low-angle sub-detectors, it is shown in [36| that it is possible to reconstruct 
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the masses of h\ and A\. In [37] it was suggested that production of A\ in association 
with charginos followed by the possibly dominant decay A\ — > 77 could offer good detection 
prospects for an almost purely singlet A±. An alternative probe is the decay T — > A17 at B 



factories [38] • A high-energy e + e~ linear collider would easily probe the scenario of dominant 
decay hi — > AiA% for < 2m^ via the recoil mass technique which is insensitive to the 
decay of h\. 

For Ma 1 > 2mfc one would have the dominant decay AiA% — > bbbb for which a LEP limit 
of Mh ± > 110 GeV was derived. In such a scenario the fine-tuning problem is not greatly 
ameliorated but detection prospects at the LHC are much better. In partonic level analyses 
it has been shown that a signal with high significance and full Higgs mass reconstruction can 



be obtained from the process pp —> Wh\ — > WA\A\ — > Wbbbb [X 31 . |35| . The main challenge 
in reconstructing the full decay chain is to retain an adequate tagging efficiency of b's in the 
low pt region where signal events are located, as shown in [35| . 

In many of the studies which are concerned with establishing a no-lose theorem the 
charged Higgs mass is taken to be very heavy M H ± > 400 GeV (e.g. the benchmark points 
in [31]). It has been known for some time that a moderately light M H ± < m t is possible in 
the NMSSM. A first detailed study appeared in [fij, and this possibility has recently been 
emphasized in 3^j. However, such a would contribute sizably to the rare decay b — > 57 



whose branching ratio has been measured and is consistent with the SM expectation. In 
the context of the NMSSM a contribution to b — > S7 from another New Physics particle 
(usually the li ghte st chargino, xt) is needed to partially cancel the large H ± contribution 



for M H ± < m t [40]. If flavour violation induced by gluinos (g) is considered [4l|, the NMSSM 
parameter space for a light H ± can be enlarged while keeping the branching ratio for b — > 57 
consistent with the measured value. This merely requires a suitable cancellation among the 
contributions from H ± , xt an d 9, the latter being of essentially arbitrary magnitude. In 
light of this possibility we do not impose the b — > 57 constraint in our numerical analysis. 

Another potentially important constraint on the scenario of Mjj± < m t comes from 
the measurement of the decay — > t^v 42|. This decay is mediated at tree-level [43[ 
by and its contribution cannot be canceled by any other new particle in the model. 
Current data excludes two re gion s in the parameter space of [Mh± , tan (3] . However, the 
non-holomorphic contribution 44| would shift the location of these two regions and thus we 
do not impose such a constraint in our analysis. Importantly, for tan f3 < 20 of most interest 
to us M H ± < m t is almost always allowed. Moreover, a recent analysis j45| shows that a 
light charged Higgs boson in the NMSSM is compatible with the constraints from b 57, 
AM q , B ± — > r ± i/ and B s — > fi + fi~ even without invoking extra sources of flavour violation 
from gluinos. 



V. CHARGED HIGGS DECAY H± -> SW AND THE PRODUCTION MECHA- 
NISM pp -> H±S, S = Ax,h t 

A. Charged Higgs decay modes — > SW 

The decay — > AW, where A is a CP-odd Higgs boson, may be sizeable in a variety 
of models with a non-minimal Higgs sector such as Two Higgs doublet models (Type I 



and II) [46|, |47|, |48| and in SUSY models with Higgs triplets [49]. Two LEP collaborations 
(OPAL and DELPHI) performed a search for a charged Higg s decaying to AW* (assuming 
rriA > 2mb) and derived limits on the charged Higgs mass [50| comparable to those obtained 
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from the search for H ± — > cs, rv. In the MSSM the decay width for if ± — > AW is very 
suppressed in most of the parameter space [H, 0] because the charged Higgs and the CP-odd 
Higgs are close to mass degeneracy 

The importance of the decays H± -> A X W and H ± -> h{W in the NMSSM was first 
pointed out in |§]. Their branching ratios may be close to 100% which can provide a clear 
signal at the LHC. Simulations of the process pp —> it followed b y t — > and H ± — > A±W 



have been performed for the NMSSM [7[, CP conserving MSSM [51| and CP violating MSSM 
13] . The partial width is given by: 



^ - A ™ = 16 ^/ j (^X, *S0 (!3) 

where A(x, y, z) = x 2 + y 2 + z 2 — 2{xy + xz + yz) is the two-body phase space function. The 
decay width of H ± — > h{W can be obtained from eq. ( {TBI by replacing cos6U by g^w^hx 
and M Al by M ftl . 

As can be seen from ( JT3|) . the decay width of if ± — > AiVK is directly proportional to cos ^ 
which is the doublet component of A\. This decay width can be substantially enhanced if 
A\ is predominantly composed of doublet fields. However, even with small doublet (large 
singlet) component of A\ it is possible that H ± — > A\W is the dominant decay mode. We 
perform a scan of the parameter space using the code [s| (NMSSM- Tools incorporates the 
LEP2 bounds for 4b and 6b final states) in order to quantify the importance of H ± — > A{W 
and #± -> h x W. 

Hereafter we assume that all scalar superparticles share the same soft mass term Msusy-, 
and the ratios of gaugino masses satisfy M\ : M 2 : M 3 = 1:2:6; the trilinear couplings are 
related to Msusy but the sign is not fixed, i.e. A tj b = ±2Msusy- We scan the parameter 
space of the model by varying the free parameters within the following region: 

A = [0,1], re =[-1,1], tan (3 = [0.2,60] , /i = [-1, l]TeV , 
A x = [-1.0, 1.0]TeV, A K = [-1.0, 1.0]TeV, 

Msusy = [0.2, 3]TeV , M x = [0.07, 3]TeV . (14) 

While varying these parameters, we take into account the experimental constraints on the 
MSSM spectrum e.g., charged Higgs mass > 80 GeV, chargino and scalar fermions > 100 
GeV. We also apply the full set of LEP constraints obtained from searches for neutral Higgs 
bosons decaying to final states like Z2b, Z4b, 6b, 6r, Z2b2r, Z4t, 262t. 

In Fig. ([T]) we display the branching ratios of W^Ai , rv and top-bottom modes. Before 
the opening of the H ± — > tb channel, the full dominance of W ± A X over rv requires light 
Ma x ^ 100 GeV, large doublet component of A\ and tan/3 not too large. Note that at large 
tan/3 pa 15 — 25, the W A\ and rv channels become comparable in size. Once the decay 
tb is open, it competes strongly with W A\ for tan/3 < 15. As can be seen from 
Fig. ([1]) upper left, the branching ratio of if ± — > W A\ is less than 90%. It is interesting to 
see also that for cos 2 9a ^ 0.05 there is not a single point with Br(H ± — > W ± A 1 ) > 50%. 
Note also that at large tan/3 > 25, it is hard for H ± — > W ± A 1 to compete with tu and 
top-bottom modes. 

The case of the analogous decays if ± — > W ± hi^ are displayed in Fig. ([2]) as a function of 
Mh± and tan/?. One can see from the upper right panel of Fig. (T5]) that W ± h\ dominates 
over rv only for moderate tan/3 < 5 and before the opening of — > tb decay, which 
strongly competes with H ± — > W ± h\ mode. 
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FIG. 1: Comparison of the branching ratios of — » {W^Ai, tu, tb} as a function of M H ± (upper 
left), cosOa (upper right), Ma 1 (lower left) and tan/3 (lower right). In all panels only points with 
Br (if* -> W ± Ai) > 50% are selected. 



^From the lower panel of Fig. (j2J) one can see that the branching ratio for —*■ W ± fi2 
can only be larger than 20% for charged Higgs mass larger than about 220 GeV. This is 
mainly due to the fact that is most of the time larger than 140 GeV. It is clear that 
both H — > W ± li2 and H ± — > tb are of comparable size except in the case of large tan (3 
where — > tb mode dominates. 

Importantly, we note that if Sf 3 ~ 1 the second sum rule in Eq. (jSJ) requires gy Vhl ~ and 
9w±uTh x ~ 0- In this case, Sf 3 ~ 1, both modes H ± — * W ± h\ )2 are suppressed and hence 
the full dominance of W ± hi^ requires small S^. 

In our numerical analysis we have explicitly checked that if h\ is predominantly singlet, 
i.e., S\ 3 > 0.9, both couplings Qvvh^ 9w ± HTh- L ~ ^.1, in accordance with this sum rule. The 
larger S\ 3 is, the smaller are the couplings g\y hx and g^± H ^ h ■ When S\ 3 > 0.9, h\ is almost 
purely singlet and even a very light hi can be allowed by LEP experimental constraints. In 
this case, both the vertices of ZZh\ and ffh\ are suppressed, as shown in Fig. ((3^). In 
this case, the hi will be the Standard Model like Higgs boson and the coupling gvvh 2 can 
be large, as indicated by the first sum rule in Eq. (|S]) and demonstrated in Fig. 

In the converse case when Sf 3 ~ — > 0.1, g VVhl and g^ v ± H ^ hx have to share the quantity 
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FIG. 2: Comparison of the branching ratios of — > {W^h\,Tv} (upper plots) and — > 
{W ziz h2,tb} (lower plots) as a function of M H ± (left) and tan (3 (right). In all panels only points 
with Br(F ± -> W ± h i ) > 20% are selected (i = 1 or 2). 

1 — Sf 3 . Since hi is dominantly doublet the coupling g w ± H ^hi can be maximal, and hence 
the branching ratio of H ± — > W ± hi can be large. 



B. The cross-sections for pp — > H^hi, pp — > W h\ and — ► ff ± A 1 in the NMSSM 

Searches for Higgs bosons at the LHC suffer from large QCD backgrounds. However, 
detailed studies have shown that multiple signals for the MSSM Higgs bosons are possible 
in a sizeable region of the plane [tan/5, M H ±] j53[. Much of these studies for the MSSM can 
be applied to the NMSSM with some caveats which were discussed in Section IV. The most 
problematic region for H ± discovery in the MSSM is for moderate values of tan/?, since the 
production mechanisms which rely on a large bottom quark or top quark Yukawa coupling 
(e.g. gb — > H^t) are least effective. Hence alternative mechanisms which could offer good 
detection prospects for at moderate values of tan j3 are desirable. 

The cross sections for the pair production mechanisms pp — > H ± Ai and pp —>■ H ± hi fall 
quickly with increasing scalar masses but for relatively light masses (< 200 GeV) they can 
provide promising signal rates which might enable their detection at the LHC. One common 
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FIG. 3: Left panel: the comparison ofm^ andm^ with respect to Sf^. Right panel: the comparison 
of the cross section of the processes pp — > W h\ and pp — ► hi. For a(pp — > H^Hi), we sum 
over a(pp — > H + Hi) and o~(pp — > H~ Hi), while for o~(pp — ► W^/ii) we sum over <r(pp — > 
and o"(pp — > For i/ie sa/ce o/ comparison, we deliberately divide o~(pp — > W hi) by three 

due to the three helicity states of massive vector boson W . 



feature is that the produced scalars enjoy large transverse momenta, which are crucial for 
the trigger and event selection. The cross section for pp — > — *> H^A was first studied 
at both the LHC and Tevatron in the CP conserving MSSM for M A > 100 GeV. The 
analogous process pp — > H ± hi for a very light hi with unsuppressed coupling hiH ± W Zf was 
studied in the 2HDM and the CP violating MSSM in [ll[ at the Tevatron. In 1(| it was 



shown that pp — » H ± h°, H ± A° can be important in specific regions of parameter space (i.e., 
very light h°, A ) in the CP conserving MSSM. 

In the NMSSM, if the coupling H ± W Zf Ai is sizeable, so will be the cross section for pp — > 
W ± —> H^Ai provided that and Ai are not too heavy. The production mechanism pp — > 
H^Ai followed by the decay H ± — > W ± Ai would give rise to a signal W-^AiAi — > Wbbbb 
or W /± y4iv4i — > Wtttt. The signature iy ± AiAi —> Wbbbb was simulated at the LHC 



in [12|] in the context of the CP violating MSSM with the conclusion that a sizeable signal 
essentially free of background could be obtained. We use NMSSM-TOOLS1.1.1 to calculate 
the mass spectrum and couplings of the NMSSM Higgs bosons, and we link CTQ6.1M PDF 
distribution to this code in order to calculate the cross sections of pp — > H Ai, pp — > H^hi 
and pp — > W^hi. All cross sections are evaluated at a scale which is the sum of the masses in 
final states and do not include next-to-leading order QCD enhancement factors (K factors) 
of around 1.2 -> 1.3 0,0. 

For our numerical analysis, we have done a systematic scan with NMSSM-TOOLS1.1.1 
Firstly, we explore the phenomenological implication of the sum rule Eq. (jSj) with Fig. 
((3b). There are several comments in order: 

1) All points in Fig. ([3b) respect the following constraint M H ± > M w , this leads to a 
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smaller cross section for a(pp — > H ± h\). 

2) As Sf 3 increases both processes are suppressed due to the decrease of the couplings 
W^W ± h 1 and W^H^hx- 

3) When Sf 3 > 0.8, hi is dominated by singlet component, therefore it can be very light see 
Fig. (j3h)- In this cases, according to sum rule Eq. ([8]), the vertex VVh\ suffers a severe 
suppression. However, some points with large a(pp — > W ± h x ) arise due to the fact that a 
very light hi is allowed. 

In Fig. (jl^i) we study the cross section of pp — > H hi at the LHC and select points which 
simultaneously satisfy the following conditions: 

a(pp -> H+hx) > 0.1 pb and Br(H ± -> W ± A 1 ) > 0.5 . (15) 

We require points in parameter space with cross sections larger than 0.1 pb as a conservative 
threshold of observability for this channel at the LHC. From the figure it is clear that 
cr(pp — * H ± hi) < 0.5 pb when the charged Higgs boson decays dominantly to W ± Ai. 

In Fig. (j4b), we study the cross section of pp — > W ± hi at LHC and select points which 
satisfy the following conditions: 

a(pp -> W ± hi) > 0.1 pb and Br(hi -> Ai Ai) > 0.5 . (16) 

The typical cross section for a(pp — > W ± hi) is around a few pb, which is considerably larger 
than a(pp — > H hi). The larger cross sections correspond to the larger branching ratios 
for hi — > AiAi (> 90%). The numerical results in Fig. (jib) are in good agreement with 



analogous results presented in (32 



In Fig. (JSj) we analyze the components of Hi and Ai. Points which satisfy the following 
condition are selected: 

Br(hi -> AiAi) > 0.5 . (17) 

As expected, when both and Ai are dominantly composed of doublet fields the region of 
light Higgs bosons is ruled out from searches for e + e~ — > Zhi — > Z2Ai — > Z4b, and 
should be heavier than around 100 ~ 110 GeV. When Ma x ^ 2mft, can be lighter than 
100 GeV due to the fact that the LEP2 sensitivity to the channel e + e~ — > ZAt was less 
robust than that for e + e~ — > ZAb. Interestingly, when both hi and Ai are mainly singlet 
and hence the vertex of VVhi is greatly suppressed, much lighter values for (< 80 GeV) 
are still allowed, as shown by points with red stars in Fig. (j5^i) and blue crosses in Fig. (j5b). 

This process pp — > H ± Ai — > H/ ± AiAi leads to the same signature as the process pp — > 
W/ii — > WAiAi — > Wbbbb. The latter has been simulated in |l3j and also offers very good 
detection prospects. We will compare the magnitude of these two distinct mechanisms which 
lead to the same Wbbbb signature. In addition, the mechanism pp — > H^hi followed by the 
decay —>■ W ± Ai would also lead to the same final state W ± A 1 h 1 — > Wbbbb. We will 
concentrate on the scenario where hi — > AiAi is large and thus hi — > bb will be kinematically 
suppressed. We will discuss the magnitude of pp H ± hi —> W ± Aihi — > Wbbbb later. 

In Fig. ([6^) we study the process pp — > H ± Ai by choosing points which satisfy the 
following conditions: 

a(pp -> H±Ai) > 0.1 pb and Br(H ± -> W ± A 1 ) > 0.5 . (18) 

It is apparent that the magnitude of a(pp — > H ± Ai) can reach a few pb and thus is within 
the detection capability of the LHC. The analysis of [l2| (for the CP violating MSSM) 
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o(pp^H ± H 1 )(PB) c<pp->W ± H 1 )(PB) 

FIG. 4: Left panel: points selected with the condition given in Eq. $15\) . Right panel: points 
selected with the condition given in Eq. I116\) . For o~(pp — > H^h\) we sum over a(pp — > H + h\) 
and cr(pp — > H~h\); for a(pp — > W^/ii) we sum over a(pp — ► W + h\) and cr(pp — > W~/ii). VKe 
s/tow the two decay modes of A\: A\ — » 65, and Ai — ► rf, which corresponds to two mass regions: 
2mt, < Ma 1 < and 2m T < Ma x < 2m?,, respectively. 




IAD 



m H . (GeV) m... (GeV) 



FIG. 5: Parameter space satisfying Br(H\ — > ^4i^4i) > 0.5 in £/ie plane [m/ ll ,m J 4 1 ]. T/ie compo- 
nents of both hi and A\ are displayed. 
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FIG. 6: Left panel: Points in the plane [o~(pp — > H^Aij^Br^^ — ► .AiW] which satisfy the 
conditions given in Eq. \18\). For a(pp — > H^A\) we sum over cr(pp — > H + A\) and o~(pp — > 
H~A\). We show the two decay modes of A\: A\ — > bb and A\ — > rf, which corresponds to 
two mass regions: 2mi, < Ma 1 < and 2m T < Ma 1 < 2m^, respectively. Right panel: the 

dependence of o~(pp — > H^Ai) on both tan/3 and cosOa are displayed. 

suggests that a(pp — > H ± Ai) > 0.1 pb with a large i?r(if ± — > W /± t4 1 ) would be sufficient 
for an observable Wbbbb signal at the LHC Most strikingly, the cross section of the process 
pp — > H ± Ai can be comparable to that of pp — > W ± hi. 

The majority of the points in Fig. ([6h) correspond to the parameter space where tan/3 
is located in the range 0.2 < tan/3 < 20. As seen in the previous section, when tan/? > 20 
the decay channel H ± — > t^Uj. (or — > tb ) will dominate over H ± — > H /± A 1 . It is 
evident from Fig. ([6^) that there are plenty of points with a{pp — > H^Ai) > 0.1 pb and 
Sr(iJ ± -> W ± A l ) > 90%. 

In Fig. ([6b) we show the dependence of a(pp — > H^Ai) on tan/3 and cosOa- The figure 
clearly shows that when A\ is mainly doublet the cross section a(pp — > H^Ai) can reach a 
few pb. Importantly, the cross section can be sizeable in the whole region 1 < tanfi < 30, and 
thus this mechanism can be applied to the region 5 < tan/3 < 20 for which if ± discovery in 
the conventional production mechanisms (which utilize the t and b quark Yukawa couplings) 
are least effective. Thus H production via pp — > H ± Ai might offer the best prospects for 
the detection of a light NMSSM charged Higgs boson in the region of intermediate tan /3. 
It is clear from Fig. (JBb) that there are no points at all with < cos 2 9a ^ 0.4, the reason 
being that such points do not satisfy the requirement a{pp — > H ± A\) > 0.1 pb. 

Fig. ([7]) shows the dependence of a{pp — > H ± hi{Ai)) on m^m^J +m H ±. Points in Fig. 
([7^) satisfy the conditions given in Eq. f|T5|) . while points in Fig. ([7b) satisfy the conditions 
given in Eq. f[T8|) . Clearly the points with large cross section correspond to the region in 
the parameter space where both and h\(Ai) are light and the couplings W Zf H ± hi and 
W T H ± A 1 are near maximal. In Fig. ([7b), it is evident that the cross section ior pp — > H ± Ai 
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FIG. 7: Left panel: the cross section of pp — > H^h\ against + M H ±; all points satisfy the 
condition in Eq. H5\). Right panel: the cross section of pp — » H^A\ against Ma 1 + M H ±; all 
points satisfy the condition in Eq. H8\). We show the two decay modes of A±: Ai — > bb and 
A\ — » rf, which corresponds to two mass regions: 2m& < M\ x < and 2m T < < 2m;,, 

respectively. 

can reach a few pb when A\ is as light as 10 GeV. In contrast, in Fig. (JT^l) one can see that 
that there are only points for m/ ll + m#± > 170 GeV which corresponds to m#± > 80 GeV 
and ft 90 GeV. The lack of sample points with large cross section for pp — > H h\ is 
due to difficulties in finding points with relatively light hi and H^ 1 (i.e., + i71h± ^ 170 
GeV) which can satisfy the experimental constraints. 

As emphasized earlier, the processes pp — > H Ax and pp — > Vh\ could lead to the same 
final state, Wbbbb or Wtttt. Hence a numerical comparison of their cross sections is of 
particular interest and is shown in Fig. (jSJ), where all points satisfy the following conditions: 

a{pp -> iZ^Ai) > 0.1 p6 and a(pp -> W ± h 1 ) > 0.1 pb . (19) 

Superimposed on Fig. and Fig. ([8b) are the main decay modes of the charged Higgs 
boson and the decay neutral Higgs boson Hi respectively. We further impose the following 
conditions: 

BriH* -> W ± Ai) > 0.5 and Br(hi -> A x Ai) > 0.5 , (20) 

and the surviving points are displayed in Fig. (jU^). Importantly, there are many points 
where the two cross sections are of comparable size. We note that for these points in Fig. 
flSt) the pseudoscalar Ai can be both R-axion like or a mixture of the three allowed basic 
axions. If the magnitude of the cross sections of both pp H A\ and pp — * Vhi are similar 
then the interference of the two channels (i.e., the same Wbbbb signature arising from distinct 
production mechanisms) should be taken into account. We have neglected such effects in 
the present study. 
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□(pp-sH* A,)(PB) ujpp^H 1 A,»(PB) 

FIG. 8: Left panel: comparison of a(pp — > H^Ai) and o~(pp — > W^hi) with two decay modes. 
Right panel: comparison of a(pp — > ff^i) and a(pp — > W^/ii) w«£/i iwo /ii decay modes. The 
dotted line corresponds to cr(pp — » W /ti) = cr(pp — > ff^^Ai). 

We now discuss whether the Wbbbb signatures can be distinguished experimentally by 
comparing the strategies adopted in [H] (for pp — > f/^yl ) and [l3[ (for — > W ± h\). 
In order to reconstruct the peak of the CP-even Higgs h -i , o ne can select events with a 
charged lepton and four tagged b quark jets as shown in [l3j. This enables both a clean 
Higgs signal with high significance and a measurement of given by the invariant mass 
of the four b quark jets, m^. The process pp — > H ± A 1 might be an irreducible background 
but presumably could be significantly suppressed with the aforementioned cut on e.g., 
m hl — 15GeV < m ib < m hl + 15GeV. 

Regarding detection of pp — > H ± A°, it was demonstrated in [l2| (for the analogous 
process pp — > H ± H± — > WH\H\ in the CP violating MSSM) that the mass of H ± can be 
reconstructed. This is achieved by defining a tranverse mass (Mt) which is a function of the 
momenta of the two secondary b jets (i.e., those originating from the decay H ± — > A{W — > 
Wbb) and the momenta of the lepton and missing energy coming from the W boson. It was 
shown that Mt is sensitive to the underlying charged Higgs mass and thus can be used for 
the determination of Mh±- The pair of b jets from pp — > W hi might be an irreducible 
background but presumably could be suppressed with a cut on Mt 

To reconstruct the peak of the light CP-odd neutral Higgs A\ one can require events with 
four tagged b jets, construct the three possible double pairings of bb invariant masses, and 
then select the pairing giving the least difference between the two bb invariant masses values 
12|. W4b signatures from the process pp —>■ W ± hi also contribute constructively to the 



reconstruction of A\. Thus we conclude that it is promising to reconstruct the peaks of the 
CP-even neutral Higgs (hi), charged Higgs (H^) and CP-odd neutral Higgs (Ai) and thus 
experimentally distinguish the Wbbbb signatures arising from the two distinct production 
mechanisms. We defer a detailed simulation to a future work. 
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FIG. 9: Left panel: comparison of a(pp — > H^Ai) and o~{pp — ► W^h\) with different A\ decay 
modes. Points are selected with the condition given in Eqs. \19\2Q\i . Right panel: comparison of 
a(pp — ► H A±) and o~(pp — > H^hi) with the same set of points. The dotted line corresponds to 
a(pp — > il Ai) = — ► H^hi); the dashed line corresponds to o~(pp — > iT Aj.) = 3a (pp — > 
H^hi); the solid line corresponds to a(pp — > H^Ai) = 10cr(pp — > H^hi). 

Finally, we also compare the cross sections of pp —>■ H ± Ai and pp — > H ± hi in Fig. ([9b). 
The points are from the same data sample used in Fig. ([9^). It is clear that a(pp — > H^Ai) is 
around one order of magnitude larger than a(pp — > H ± h 1 ), and the underlying reason is that 
Mh x > 2Ma 1 - Consequently, pp — > H hi — > W ± Aih\ — > Wbbbb will also be suppressed and 
can be safely neglected. Another interesting feature from Fig. ([Hb) is that points satisfying 
the conditions listed in Eq. (I2UI) lead to A\ composed mainly of the doublet fields. The 
conditions in Eq. (|20|) together with the dominance of h\ by doublet component (small Sis) 
can give large cross sections for both channels. 



VI. CONCLUSION 

In summary, we have studied the phenomenology of light charged Higgs bosons in the 
framework of NMSSM. We performed a comprehensive study of the magnitude of the branch- 
ing ratios for the decays H ± — > W ± Ai and — > W ± h\ (first considered in (6[). It was 
shown that such decays can dominate over the standard decays — > t v and E ± tb 
both below and above the top-bottom threshold. This is due to the fact that A\ can have 
a large doublet component and small mass. Large branching ratios for H ± — > W Ai and 
— > W^hi would affect the anticipated search potential for at the LHC. 

We also studied the production process pp — > H ± A 1 and showed that sizeable cross 
sections (> 1 pb) are possible. We compared the magnitude of the cross sections for both 
pp — > H A\ and the Higgsstrahlung process pp — > W hi and showed that they can be 
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of similar size. If and h\ decay via H^ 1 — > W ± A\ and hi — > A\A\ respectively, the 
above two processes would lead to the same final state, Wbbbb or Wtttt. We stressed 
that the interference term for Wbbbb and Wtttt might not be negligible and should be 
taken into account in any simulation study In particular, the signature Wbbbb affords 
promising detection prospects at the LHC and we discussed how to distinguish the distinct 
contributions from pp — * H ± A 1 and pp —> W ± h\ by using appropriate cuts. 

It is known that intermediate values of tan/3 (e.g., 5 < tan/5 < 20) are most problematic 
for discovery of H ± at the LHC since the H^tb Yukawa coupling (which is employed in 
the conventional production processes) takes its lowest values. In such a region the process 
pp — > H ± A\ can have a sizeable cross section if m#± + rriA 1 < 200 GeV. Therefore we 
propose pp — > H ± Ai as a unique mechanism to probe the parameter space of intermediate 
tan (3 and light charged Higgs boson in the NMSSM. 
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